Steady Stokes flow with long-range correlations, fractal Fourier spectrum, 

and anomalous transport 
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We consider viscous two-dimensional steady flows of incompressible fluids past doubly periodic 
arrays of solid obstacles. In a class of such flows, the autocorrelations for the Lagrangian observables 
decay in accordance with the power law, and the Fourier spectrum is neither discrete nor absolutely 
continuous. We demonstrate that spreading of the droplet of tracers in such flows is anomalously 
fast. Since the flow is equivalent to the integrable Hamiltonian system with 1 degree of freedom, this 
provides an example of integrable dynamics with long-range correlations, fractal power spectrum, 
and anomalous transport properties. 
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The last two decades have witnessed growing interest 
in theoretical and applied aspects of chaotic advection 
in steady or temporally periodic velocity fields of hydro- 
dynamical flows 0, 0]. The contrast between the triv- 
ial behavior of the field and the complicated dynamics 
of individual fluid particle is due to the fact that the 
field is usually described by means of Eulerian variables 
in the laboratory reference frame, whereas for the fluid 
particle the Lagrangian description is more appropriate, 
since its reference frame is carried by the flow. An Eule- 
rian observer performs measurements in the fixed point 
of the physical space and registers a steady or periodic 
process. For a Lagrangian observer, physical space turns 
into phase space, streamlines become phase trajectories, 
a "measuring device" is advected along these trajectories 
and explores different regions of the flow pattern. There- 
fore the motion of the fluid particle can be chaotic in 
spite of the time-independent state at any given place. 
A consequence of this dynamics is intensive stirring and 
mixing. Examples of chaotic advection are widespread, 
from microscopical flows of nanotcchnology and cellular 
biology to large-scale phenomena in astrophysics andgeo- 
physics; for a list of references see a recent review yj- 

To exhibit Lagrangian chaos, a flow of incompress- 
ible fluid must be three dimensional, time-dependent, or 
both. In two-dimensional setup, there is no place for ex- 
ponential divergence of streamlines, and in conventional 
examples the motion of a typical fluid particle is periodic 
and well correlated. Below, we discuss a class of steady 
two-dimensional flows for which the correlations decay 
and power spectra are fractal. 

Consider a slow steady flow of a viscous incompress- 
ible fluid past the square lattice of parallel solid circular 



cylinders. We restrict ourselves to the "transverse" case 
when the component of the flow along the cylinder axes 
is absent. This geometry is typical, e.g., for the cooling 
flows past the rods of the reactor; on the other hand, it 
serves as the starting point for the analysis of fluid mo- 
tions in composite and porous media. Usual assumptions 
of the creeping flow allow one to neglect nonlinear terms 
in the Navier-Stokes equations and reduce them to the 
Stokes equations for the velocity V and pressure p: 



-Vp + rjAV = 0, div = 



(1) 



where 77 is the viscosity of the fluid. Boundary conditions 
enforce periodicity on the borders of the elemental cell of 
the lattice, as well as vanishing velocity on the border dT 
of the cylinder T. If the axes x and y are aligned along 
the axes of the lattice, and the cell size is chosen as the 
length unit, the boundary conditions are: 

V(0,y)=V(l,y), V(M)=V(x,l), V| ar = (2) 

Formulation of the boundary problem is completed by 
prescribing the components of mean flow along both axes: 



V x dy = a, 



Vydx 







(3) 
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Equations Q-0 were treated in terms of periodic 
Green's function by H. Hasimoto who evaluated the drag 
on a cylinder |4j ; the fundamental solution was expressed 
in terms of series containing elliptic functions [j| . Perme- 
ability of the array of cylinders was estimated in 0]. As 
every two-dimensional flow of incompressible fluid, the 
problem can be cast into a Hamiltonian frame by intro- 
ducing the stream function ^(x,y): V x = d"i>/dy and 
V y = —d^/dx. Motions along the streamlines are gov- 
erned by the integrable Hamiltonian system with 1 degree 
of freedom in which V x and V y are canonical variables, 



2 



and *B plays the role of energy. In terms of W Eq. (JTJ 
turns into the biharmonic equation AA^P = 0. No-slip 
boundary conditions require d^f /dx\ Qr — /dy\ dT = 0, 
whereas conditions for mean flow and periodicity are en- 
sured by putting ^(x, y) — ay — 0x + &(x,y) where $ 
has period 1 in both arguments. 

We solved the biharmonic equation numerically, im- 
posing periodic boundary conditions upon a polar grid 
with 600x1200 nodes; for tracer dynamics, the 2nd order 
interpolation scheme was used. Figure ^a) presents the 
flow pattern with inclination a/[3 = (y/B — l)/2 and ra- 
dius of the cylinder 1/6. 



(b) 
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FIG. 1: (a) Flow pattern with a/0 = (V§ - l)/2; (b) 
divergence of passage time for streamlines near the cylin- 
der. Crosses: numerical values, straight line: power law 
~\x-x \- 1/2 



Since the flow is steady, in every point of the laboratory 
reference frame the values of velocity and pressure are 
time independent. The picture looks different in a frame 
attached to a neutral tracer carried by the flow. We 
characterize it in terms of normalized autocorrelation, 
defined a time-dependent observable as 



C(r) 



+ t))- TO' 



(4) 



where averaging over t is performed. Since autocorrela- 
tion is the Fourier transform of the power spectrum (the 
Wiener-Khinchin theorem), its properties are related to 
the nature of the Fourier spectral measure. For peri- 
odic and quasiperiodic dynamics, C(r) displays series of 
peaks with unit or asymptotically unit amplitude at arbi- 
trarily large values of r; this corresponds to the discrete 
(pure-point) power spectrum. In case of chaotic behav- 
ior (impossible in two-dimensional dynamics) autocorre- 
lation decays exponentially, and the power spectrum is 
absolutely continuous with respect to the Lebesgue mea- 
sure. 

Qualitative features of the autocorrelation are indepen- 
dent of the choice of a generic observable. We take the ab- 
solute value \ V\ of the velocity and measure it for the fluid 
particle moving along the streamline of the flow. The plot 
of C(t) is presented in Fig.Efa). The largest plotted val- 
ues of r correspond to the passage through, on the aver- 
age, 3 x 10 4 lattice cells in the y direction. We observe 



that autocorrelation decays, in spite of the ordered lam- 
inar structure of the flow. The decay is slower than ex- 
ponential; as seen in Fig.^b), the highest peaks of C(r) 
which form a log-periodic lattice, decrease in accordance 
with the power law C(r) ~ t~ - 28 . Power-law decay of 
autocorrelation ("long-range correlations") is abundant 
in physics of critical and nearly critical states and has 
been recovered in many natural phenomena, from DNA 
sequences 0, 0, U3 to atmospheric variability [Hj . How- 
ever, to our knowledge, such decay has been reported 
neither in the context of steady two-dimensional flows, 
nor in a more general context of integrable autonomous 
Hamiltonian systems with 1 degree of freedom, where 
the correlations are usually not supposed to decay at 
all. Our results show that integrability does not guar- 
antee against the decay of correlations. Properties of au- 
tocorrelation characterize both qualitatively and quan- 
titatively the set which supports the Fourier spectral 
measure. A useful tool is the integrated autocorrelation 
Cint(i) = t 1 Jo C{t) 2 <1t: if C- m t(t) vanishes in the limit 
t — > oo, power spectrum does not include the discrete 
component. Furthermore, C- m t(t) ~ t~° 2 where Dp. is 
the correlation dimension of the spectral measure [12j. 
The decay of C- m t(t) over several orders of magnitude in 
Fig-GI c ) indicates the absence of the discrete component 
in the Fourier spectrum of the tracer velocity. Estima- 
tion of the curve slope yields D2 ~ 0.82; hence the set 
which supports spectral measure, is fractal. Apparently, 
in this intermediate situation between chaos and order 
the Fourier spectrum of velocity is singular continuous: 
neither discrete, nor absolutely continuous. 

Origin of this unconventional state can be understood 
if the viscous flow is viewed as a dynamical system. 
From this point of view, Eqs Q - © describe the area- 
preserving vector field x = V x (x,y), y = V y (x,y) on a 
2-torus whose rotation number a equals a/0. Peculiar- 
ity of this field is its identical vanishing on the cylinder 



1 

1° 
-1 



. (a) 



1 10 10 10 10 



o 



10-' 
10- 2 



S 10"' 
10-" 



(b) 


1 1 1 x ° 28 


1 10 10 2 10 3 10 4 

X 






- (c) 


~~~~~ -^^-^t^ 82 ^^ 


1 1 10 10 2 10 3 10 4 



§10 
10 



(d) 



o 



110" 

10 J 



10 


10 2 


10 3 10 4 


-Mil 






Jllllii 


ill 




1 10 


10 2 

T 


10 3 10 4 


- (f> 




t -0.9 


1 10 


10 2 


10 3 10 4 



FIG. 2: Autocorrelation and integrated autocorrelation, (a)- 
(c): Stokes flow QJl-Q; (d)-(f): special flow ©. Dashed 
lines: power laws. 
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border dT: no-slip boundary conditions turn each point 
of dT into the fixed point. Below, we deal with irrational 
values of a. Most of the existing work on dynamics on 
2-tori, starting with the classical paper in [13j, refers to 
fields without fixed points. Consider a mapping of a line 
transversal to the flow (e.g. the boundary y — of the 
square) onto itself. In absence of fixed points the "return 
time" T re t(x) needed for one iteration of this mapping 
(a turn around the torus) is bounded. As a consequence, 
dynamics of the flow is adequately reproduced by the cir- 
cle map x — > (x + <j) mod 1. For a passive tracer in such 
quasiperiodic flow, the Fourier spectrum is discrete, and 
the autocorrelation does not decay. 

Presence of stagnation points in which the velocity van- 
ishes qualitatively changes dynamics. In this situation 
T re t (a;) as a function of x diverges, and equivalence with 
the Poincare map is invalidated: the flow with contin- 
uous power spectrum can induce the map with purely 
discrete spectrum 0. Near the streamline which leads 
into a generic isolated stagnation point, the return time 
has a logarithmic singularity. Such stagnation points ex- 
ist in forced viscous flows with doubly periodic arrays of 
steady vortices, and the Fourier spectra in these flows are 
fractal 0. Logarithmic singularities of return times are 
too weak to make the flow mixing . Stronger, power- 
like singularities can produce mixing |l7j : however, for 
isolated stagnation points such singularities occur only 
in degenerate, structurally unstable cases. 

In the flow pattern of the whole curve dT can 

be viewed as a structurally stable family of stagnation 
points. Fluid particles exhibit the strong slowdown 
during their motion along the edge of the obstacle. The 
collective effect of this continuum of equilibria manifests 
itself in the power-law divergence of return time: let 
(xq, 0) lie on the streamline which separates the particle 
paths passing the obstacle "on the left" from the paths 
passing "on the right," then r rct (x) ~ 1/ y/\x — xq\ 
[Fig. nj)]. The existence of this singularity ensures the 
decay of correlations and creates continuous component 
in the power spectrum. 

Let us proceed from advection of individual passive 
tracer to the transport of the ensemble of tracers by the 
flow d)-©. Figure |21 shows the temporal evolution of 
initially compact droplet formed by 10 4 fluid particles. 
Nonisochronicity of rotations around the torus stretches 
the droplet. The strongest effect is observed during the 
passages in which the obstacle is sandwiched between 
the streamlines: some particles hover in the vicinity of 
the obstacle for the very long time, whereas the particles 
which pass at a larger distance, are long gone. Because 
of irrational rotation number, each streamline is dense on 
the torus; hence, each particle passes arbitrarily close to 
the cylinder where its velocity becomes arbitrarily small. 
As a result, after a hundred of passages the particles are 
spread over virtually the whole area of the square. Nat- 
urally, this mixing is slower and less efficient than mix- 




FIG. 3: Cloud of passive tracers carried by the flow. 

(a) t=0; (b) t=0.5; (c) t=2; (d) t=5; (e) £=20; (f) £=100. 

ing enabled by chaotic advection in the situation of La- 
grangian chaos 0,01 • Nevertheless, ultimately it leads to 
the same result . 

To quantify the rate of spreading of the droplet on 
the infinite plane, we introduce the displacement £ y (t) = 
y(t) — y(0) and estimate the time growth of the variance 
C 2 W = — (£yW) 2 by averaging over all initial 

positions within the unit cell. Figure0]shows that growth 
of £(i) is described by the power law whose exponent 
exceeds the "normal diffusion" value 1/2. 

For larger values of t the computed values of ((t) un- 
dershoot the expected power law. We assign this discrep- 
ancy to insufficient statistics: slow saturation in the de- 
pendence of £ on the length of the sample over which the 
averaging is made. As t grows, the main contribution into 
£(t) is made by rare "large events": passages extremely 
close to the obstacle. Apparently, the longest numerically 
affordable to us trajectories with 3x 10 6 turns around the 
torus (t ~ 2.8 x 10 6 ) are still insufficient. In order to ex- 
tend the range of t in which the estimates of are 
reliable, we have used the model which is computation- 
ally much less expensive: the "special flow" (flow over 
the mapping) , often employed in ergodic theory 0| • For 
the map u — » f(u) and the function t(u), this dynamical 
system is introduced on the part of the plane u, v defined 
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FIG. 4: Growth of variance in the ensemble of particles trans- 
ported by the flow. Crosses: Stokes flow Q-J3J; pluses: spe- 
cial flow JSJ with G = 0.3; straight line: power law t 0,6 . 
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by < v < t(u), as follows: the orbit starts at t = in 
the point (uq,vo), and moves with unit velocity parallel 
to w-axis until reaching the border t(uq); from there it 
instantaneously jumps to the point (f(u),0), begins the 
next stage of motion along v, and so on. Obviously, on 
the u axis f(u) is the Poincare mapping, and t(u) is the 
return time. The trajectory is described by the expres- 
sions 

n-i 

ut = / n (uo), v t = v +t-J2 T (f j (uo)) (5) 

3=0 

with integer n = n(t) uniquely defined by inequalities 

n— 1 n 

J2r{f(uo))<v + t<Y,r{fM) (6) 

j=0 j=o 

The special flow for our hydrodynamical problem, com- 
bines the circle map f(u) — (it + a//3) mod 1 with the re- 
turn time t{u) = C/ ' <J\u — 1/2|. Properties of autocor- 
relation for the flow ©, presented in Figs. EId)-Etf), are 
in qualitative and good quantitative accordance with the 
respective characteristics of the Stokes flow: autocorrela- 
tion decays as the power law, and the Fourier spectrum 
is fractal. To evaluate the transport in J^Jl we used the 
lift of f(u) onto the real axis. Estimates obtained from 
10 10 iterations of / are shown in Fig. 0] by pluses. Again, 
we see the anomalously fast growth of the variance. 

These unusual transport properties are related to the 
distribution of return times. For T large, the propor- 
tion of returns with duration between T and T + dt 
decreases as ~ T~ 3 dt; the second moment of this dis- 
tribution (mean square of the return time) is infinite. 
This is reminiscent of Levy flights: random walks with 
power-law statistics and infinite moments. Widespread 
in physics, Levy flights exhibit long-range correlations 
and anomalous transport 0, Of course, the com- 
pletely deterministic motion of a tracer carried by a flow 
(|T|) is not a random walk; the sequence of return times is 
rigidly prescribed by the rotation number. Nevertheless, 
as we see, already the presence of the power-law tail T~ 3 
appears to be sufficient for the onset of typical hallmarks 
of the Levy flights. 

In dynamics of nonintegrable Hamiltonian systems, 
anomalous diffusion is ascribed to sticking of trajectories 
near the Kolmogorov-Arnol'd-Moser surfaces and Can- 
tori [2llE^|. In fluid dynamics, Levy flights were exper- 
imentally observed in time-dependent flows with coher- 
ent structures (jets, vortices etc.), which serve as kind 
of traps for tracer particles |U 0j 113 • O ur example 
shows that similar anomalously fast transport can take 
place also in the framework of the two-dimensional time- 
independent flow. Moreover, both basic ingredients, irra- 
tional rotation number and vanishing velocity along the 
whole borderline, are generic within the considered class 
of Stokes flows. Of course, for a "random" rotation num- 



ber a/ (3 the peaks of the autocorrelation would not form 
a regular lattice. 

In conclusion, we have studied the spectral and trans- 
port properties of a plane steady viscous flow past an 
array of solid obstacles. On reducing the problem to the 
time-independent integrable Hamiltonian system with 1 
degree of freedom, we find fractal Fourier spectrum, long- 
range correlations and anomalously fast "superdiffusivc" 
transport reminiscent of Levy flights. 
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